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One of the principal models of magnetic sensing in migratory birds rests on the quantum spin- 
dynamics of transient radical pairs created photochemically in ocular cryptochrome proteins. We 
consider here the role of electron spin entanglement and coherence in determining the sensitivity of 
a radical pair-based geomagnetic compass and the origins of the directional response. It emerges 
that the anisotropy of radical pairs formed from spin-polarized molecular triplets could form the 
basis of a more sensitive compass sensor than one founded on the conventional hyperfine-anisotropy 
model. This property offers new and more flexible opportunities for the design of biologically inspired 
magnetic compass sensors. 
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The biophysics and biochemistry that allow birds to 
sense the direction of the geomagnetic field (25-65 f/T) 
are for the most part obscure. One of the two cur- 
rently popular hypotheses (the other involves biogenic 
iron-oxide nanostructures [T]) is founded on magneti- 
cally sensitive photochemical reactions in the retina [2]. 
It is thought that photo-induced radical pairs in cryp- 
tochrome, a blue-light photoreceptor protein, may con- 
stitute the primary magnetic sensor [3J 0] and a variety 
of supporting evidence has accumulated over the last few 
years (reviewed in [5 8J). If this mechanism proves to be 
correct, it will incontrovertibly come under the umbrella 
of 'quantum biology' [9], as an instance of Nature us- 
ing fundamentally quantum behaviour - in this case the 
coherent spin dynamics of radical pairs - to achieve some- 
thing that would be essentially impossible by means of 
more conventional chemistry. For this reason, the avian 
magnetic compass has attracted the attention of quan- 
tum information theorists and others wishing to under- 
stand the role played by spin-entanglement and to deter- 
mine whether the techniques of quantum control could 
shed light on this intriguing sensory mechanism [10H13] . 

A fundamental property of radical pairs that allows 
sensitivity to magnetic interactions orders of magnitude 
smaller than k^T is that their chemical transformations 
conserve electron spin. Radical pairs are therefore cre- 
ated with the same spin-multiplicity (singlet or triplet) 
as their precursors. Owing to electron-nuclear hyper- 
fine (HF) interactions, neither singlets nor triplets are, 
in general, eigenstates of the spin Hamiltonian. Con- 
sequently, the radical pair starts out in a non-stationary 
superposition which evolves coherently at frequencies de- 
termined by the HF interactions and also, crucially for 
a magnetic sensor, by the electronic Zeeman interactions 
with an external magnetic field [5]. Spin decoherence 
and spin relaxation can be slow enough to allow even an 
Earth-strength magnetic field to modulate the spin dy- 
namics and hence alter the yields of the products formed 
by spin-selective reactions. The anisotropy of the HF in- 



teractions leads to anisotropic reaction yields and hence, 
in principle, a magnetic direction sensor (14( 115) . 

The singlet state - the initial state of the radical pairs 
formed photochemically in cryptochromes [H [TS] - is en- 
tangled: 

|S)(S| = i|ai/? 2 )(«ift| + \\M{PiOL 2 \ (1) 
-\\aiP2){P\a 2 \ - ^\f3ia 2 )(a 1 f3 2 \ 

(a and f3 are the m$ = ±^ spin states of the two unpaired 
electrons). But other initial states are also known to re- 
sult in magnetically sensitive chemistry [T7] : do they too 
need to be entangled or is it sufficient if they are 'merely' 
coherent? Or is neither entanglement nor coherence nec- 
essary for a magnetic compass? 

Questions such as these have been addressed in two re- 
cent papers. Briegel and his group noted that randomly 
generated separable (i.e. not entangled) initial states 
could result in reaction product yields more anisotropic 
than those produced from an initial singlet state under 
the same conditions [10] . The other study, by Benjamin 
and colleagues, reached similar conclusions by analysing 
model radical pair systems, finding significant product 
yield anisotropies for the separable initial state [TT] 

||S)(S| + i|T )<T | = Hax&Xcu&l (2) 

in which T is the m§ — triplet spin state. 

Here we examine the role of initial entanglement 
and attempt to clarify the various sources of magnetic 
anisotropy that might form the basis of a radical pair 
compass sensor in birds. 

Initial radical pair states. We start by identifying 
chemically feasible initial electron spin states. Gemi- 
nate radical pairs are normally formed by spin-conserving 
chemical reactions so that at the moment of their creation 
they are either pure singlet, described by the initial elec- 
tron spin density matrix po = po(S) = |S)(S|, or pure 
triplet p = po(T) = | (l — |S)(S|). Occasionally, singlet 
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and triplet formation channels operate in parallel |18j , in 
which case po is a weighted sum of po(S) and /5o(T), i.e. 
of |S)(S| and 1: 

po = M)(S) + (1 - M)Po(T) (3) 
= i(4 At -l)|S)(S| + |(l- M )l 

Eq. ([3]) is also appropriate for 'F-pairs' |17j formed from 
radicals with uncorrelated spins (i.e. p = |). The op- 
erators |S)(S| and 1 and their linear combinations are 
invariant to rotations in the electron spin-space, mean- 
ing that all states that can be written in the form of 
Eq. ([3| are isotropic. Any p that cannot be so written 
is necessarily anisotropic. 

Significantly different initial states can occur when the 
radical pair comes from a molecular triplet precursor 
formed by intersystem crossing (ISC). This route is com- 
mon in photochemical reactions of the general type: 

AB S [AB]* T [AB] * rcaction > T [A* B'] (4) 

in which the final step that creates the triplet radical pair 
could be homolysis (as shown) or inter- or intramolecu- 
lar electron transfer, hydrogen atom transfer, etc. The 
formation of T [AB]* from S [AB]* requires the creation of 
spin angular momentum at the expense of orbital angular 
momentum. This process is mediated by spin-orbit cou- 
pling and is anisotropic in the molecular frame |JB . That 
is, the three triplet sub-levels of T [AB] * are differentially 
populated leading to a spin polarization in the molecular 
frame that is passed to the radical pair on its formation. 
In an appropriately chosen molecular axis system, the 
initial state of the radical pair may be written: 

p = PM( T i\ ( 5 ) 

q=x,y,z 

Anisotropic ISC is known to be responsible for a variety 
of spin-chemical and spin-polarization phenomena |17[ 

CS5H22- 

Aside from linear combinations of Eqs ([3| and ([5]), 
there are no other commonly occurring initial conditions 
for radical pairs subject to weak magnetic fields. 

Minimal radical pair model. Insights into the spin dy- 
namics of the various initial states just identified can be 
obtained from a minimal model |23j comprising two elec- 
tron spins one of which is coupled to a spin-Vb nucleus 
(e.g. 1 H). The HF interaction is either isotropic or axi- 
ally anisotropic according to the value of a dimensionless 
parameter, a [14]. Two cases are considered specifically: 
a = (isotropic) and a = — 1 (the anisotropic interac- 
tion that results in the largest reaction yield anisotropy 
for this 3-spin system (T3J). To account for the chemi- 
cal reactivity of the radical pair, we adopt the 'exponen- 
tial model' [23] in which singlet and triplet states react 
spin-selectively with the same first-order rate constant, 
k, to form distinct products. The quantum yields of 
these competing reactions are calculated using standard 



methods [TH [S3] (outlined in the Appendix). The two 
quantities of interest are $s, the fractional yield of the 
product formed via the singlet pathway, referred to here 
as the 'reaction yield', and A$g, the magnitude of its 
anisotropy: A<i>g = max {$s} — m i n {^s}- The variation 
of $s with the orientation of the radical pair in a 50 pT 
magnetic field is the basis of the compass sensor. 

To begin, we choose the isotropic initial condition in 
Eq. ([3| together with an anisotropic HF interaction (a = 
-1). In the not unrealistic limit, \a\ > ui » k [14l 124] : 

$ s = j + n(4/*-l) cos 2 9; A$ s = i|4 M -l| (6) 

where a is the isotropic HF coupling constant, uj is the 
strength of the magnetic field, and 9 is the angle between 
the symmetry axis of the HF tensor and the magnetic 
field vector. <&g is anisotropic, and therefore potentially 
suitable as a magnetic compass, except when the initial 
state is a statistical (Vi : %) mixture of singlet and triplet 
(p = |). The maximum anisotropy (A$s = \) occurs 
when the initial state is pure singlet {p = 1); for a pure 
triplet initial state (p — 0), A$s is smaller by a factor 
of three. These results were verified by exact numerical 
simulations (see Appendix). 

To quantify the entanglement of the various initial elec- 
tron spin states considered here, we use the 'concurrence' 
C(po) proposed by Wootters [25 for a two-qubit density 
operator. For the initial condition in Eq. C(po) is 
2/i— 1 when p> \ and zero when p < \ (see Appendix). 
Thus, a singlet-triplet mixture must contain more than 
50% singlet for the initial state to be entangled. The 
pure triplet state (p — 0) is not entangled, but as we 
have just seen it gives rise to a significantly anisotropic 
reaction yield. 

We now turn to a different initial condition, a linear 
combination of Eq. ([3| (with p — 0) and Eq. ^ (with 
Px=P y = 0; Vz = 1): 

p = ry|S)(S| + (l-ry)|T z )(T,| (7) 

i.e. an anisotropic mixed singlet-triplet initial state in 
which the triplet component is 100% polarized along the 
molecular z-axis. In the same limit as before (|a| »w» 
fe), but now for an isotropic HF interaction: 

$ s = §- ±(1-77) sin 2 0; Ad> s = i(l - »?) (8) 

where 9 is now the angle between the triplet polarization 
axis (z) and the magnetic field vector. The anisotropy is 
maximised when r\ = (pure |T Z ) triplet, A$g = \) and 
is at a minimum when 77 = 1 (pure singlet, A$g = 0). 
Once again, these expressions were confirmed by numer- 
ical simulations (see Appendix). We note that Eqs ^ 
and ([8]) predict identical maximum directional responses. 

The reaction yield is isotropic when 77 = 1 because then 
both the initial state |S) (S| and the spin-Hamiltonian are 
isotropic. The angle-dependence in Eq. ([8| clearly arises 
because the spin dynamics depend on the direction of the 



magnetic field with respect to the quantization (z) axis of 
the initial |T Z ) state [26]. The concurrence of the density 
operator in Eq. ^ is 2r] — 1 when r\ > | and 1 — 277 when 
77 < |. Pure singlet and pure |T Z ) triplet thus have the 
same degree of entanglement but lead to very different 
A$ s . 

Hitherto we have taken the reaction rates of the singlet 
and triplet states (kg and fcr) to be identical. Once this 
restriction is lifted, it is even possible to have magnetic 
field effects when the initial state is a statistical mixture 
of singlet and triplet: p = j/5 (S) + §p (T) = \i. To 
illustrate this point, simulations for the minimal radical 
pair with an anisotropic HF coupling are included in the 
Appendix. A$s is non-zero except when fcs = That 
is, a radical pair can exhibit magnetic compass proper- 
ties even when its initial electron spin state is neither 
entangled nor coherent. In this case the coherence arises 
during the spin evolution as a result of the differential 
reactivity of the singlet and triplet states. 

Relation between compass properties and entangle- 
ment. A complex picture emerges from these simple 
considerations. Entangled initial states can give small 
or zero reaction yield anisotropy. Non-entangled ini- 
tial states can lead to appreciable anisotropy. With two 
sources of anisotropic reaction yields - the initial state 
and the HF interactions - it is tricky to assess whether 
entanglement, or coherence in a given basis, is essen- 
tial for magnetic compass action. For example, replac- 
ing po = |S)(S| (Eq. Q) by p = §|S)<S| + ||T )(T | 
(Eq. ([2])) not only removes the initial entanglement, and 
the coherence in the {\otif32), \P1a2)} basis, it also intro- 
duces anisotropy that was not present in |S)(S|. Simi- 
larly, most randomly chosen initial states are anisotropic 
and some will give a larger A$s than does |S)(S| under 
identical conditions. In short, it appears that initial en- 
tanglement is not a particularly helpful concept when as- 
sessing the sensitivity of a radical pair compass; nor is it 
straightforwardly illuminating to consider the behaviour 
of artificial initial states. 

A radical pair compass based on initial-state 
anisotropy. The above considerations suggest an 
alternative compass design in which the directionality 
comes from the initial condition rather than the HF 
interactions. In the minimal model, the initial state 
that gives the largest reaction yield anisotropy is 
Pa = |T g )(T g | where q = x,y,z (see Appendix). We 
therefore compare |T 9 )(T g | with |S)(S| using exact 
numerical simulations (see Appendix). The possibility 
that spin-polarized triplet radical pairs might offer 
some advantage over singlets has been noted before but 
without realistic suggestions for the chemical origin of 
such initial states [26] . 

Figure [l] shows the reaction yield anisotropy of a radi- 
cal pair inspired by the flavin adenine dinucleotide radi- 
cal, FADH*, formed photochemically in cryptochromes 
|27j . One radical contains 1 H and 14 N nuclei with 
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FIG. 1. Reaction yield anisotropy, A$s, calculated (see 
Appendix) for a radical pair in which one radical contains a 
1 H nucleus (spin-Vi) and a 14 N nucleus (spin-1). k — 10 6 s _1 
and uj — 50 fiT. The HF coupling parameters (in mT) are: 
a H = -0.8; T n , xx = 0.8 8; T H , TO = -0.6 8; T H , ZZ = -0.2 8; 
a N = 0.4; T K ,xx = —0.5 8; Tn,yy = —0.5 8; T N , Z2 = 1.0 8. 
po = |S)(S| (black) and p = |T W )(T„| (green). Also shown 
are representations of the hyperfine tensors for 8 = (left) 
and 8=1 (right). 

isotropic HF couplings approximately equal to those of 
the proton and nitrogen (H5 and N5, see appendix) in 
the central ring of the tricyclic isoalloxazine ring system 
of FADH* (these being the two largest HF interactions 
in FADH* [28 ). The anisotropic components of the two 
interactions were also modelled on FADH*, but with a 
uniform scaling by a factor of 5, in the range 0.001 — 1.0. 
For the smaller values of 8, the spin-Hamiltonian is es- 
sentially isotropic. When the initial state po is a 100% 
spin-polarized triplet, A$s has significant magnitude for 
all values of 8. In contrast, when p — |S)(S|, A$s is 
essentially zero until the HF tensors become significantly 
anisotropic (8 rj 0.1). By the time the HF anisotropy is 
comparable to that in FADH* (i.e. S w 1.0), both ini- 
tial states give very similar directional responses to the 
50 pT applied magnetic field. This suggests that a spin- 
polarized triplet geminate radical pair with isotropic HF 
interactions could operate as a compass sensor just as 
well as an initial singlet state with anisotropic HF inter- 
actions. 

Indeed, there are circumstances in which, other things 
being equal, the anisotropy of the initial state might 
offer a more sensitive compass than one based on HF 
anisotropy. Biologically plausible radical pairs are likely 
to have many magnetic nuclei (mostly ^ and 14 N) with 
differently aligned HF tensors. Simulations suggest that 
the directional information potentially available from in- 
dividual HF tensors tends to be scrambled in a multi- 
nuclear radical pair, resulting in a greatly reduced A$s 
(see Appendix). A simple illustration of this effect is 
given in Fig. [2] which shows simulations of the reaction 
yield anisotropy for a spin system in which one of the rad- 
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FIG. 2. Reaction yield anisotropy, A<3>s, calculated (see 
Appendix) for a radical pair in which one radical contains 
four 1 H nuclei, all of which have axially anisotropic HF in- 
teractions with a — 0. The symmetry axes of the four HF 
tensors are directed towards the vertices of a tetrahedron. 
Three of the tensors have principal values: Tn = T22 = —1.0, 
T33 = 2.0 mT. The fourth is identical apart from a uniform 
scaling of the principal values by a factor 8. k = 10 6 s _1 and 
uj = 50 fiT. po = |S)(S| (black) and p = jT^) <T^ | (green). 
Also shown are representations of the hyperfine tensors for 
8 = 0.5, 8 = 1.0, and 8 = 1.5. 



icals contains four spin-V2 nuclei with tetrahedrally dis- 
posed axial HF tensors. When all four tensors are identi- 
cal (5 = 1), the reaction yield anisotropy for po = |S)(S| 
vanishes, by symmetry However, when the symmetry 
is reduced to C% v , by scaling the principal components 
of one of the HF tensors by a factor 5, the value of 
A$s increases but does not approach that afforded by 
Po = |Ta;)(T x | until |log 10 (5| reaches ca. 1.0. Thus it 
appears that the compass properties of a radical pair 
with many mutually cancelling HF interactions could be 
'rescued' by having a triplet, rather than a singlet, ini- 
tial condition, provided the triplet is spin-polarized by 
anisotropic intersystem crossing. 

Discussion. Having identified the initial spin-states 
in which radical pairs may be formed by chemical reac- 
tion, we revisited earlier attempts to determine the im- 
portance of entanglement and coherence as determinants 
of the anisotropic responses of radical pair magnetorecep- 
tors. It appears that the use of artificial initial spin-states 
for this purpose is somewhat confounded by their intrin- 
sic anisotropy, the effects of which may dominate the 
anisotropy conferred by the HF interactions. From these 
considerations it emerges that the anisotropy of radical 
pairs formed from spin-polarized molecular triplets could 
form the basis for a magnetic compass that is more sen- 
sitive than one based on the conventional HF-anisotropy 
model [2] in particular when the HF couplings are not 
strongly anisotropic or when the individual effects of mul- 
tiple HF anisotropies tend to counteract one another. 

Would a triplet radical pair compass be compatible 
with cryptochrome as the primary magnetoreceptor? In 
the cryptochromes investigated hitherto (bacterial p?9"] . 



plant [3] and frog [IB]), flavin-tryptophan radical pairs 
are formed as singlets. However, avian cryptochromes 
may behave differently, and there are precedents for 
triplet radical pairs in other flavoproteins [3TJI |3"T] . Su- 
perficially, it appears that flavins may be suitable for an 
initial triplet-state compass: intersystem crossing in both 
flavin mononucleotide and riboflavin at near-neutral pH 
results in fractional populations of the zero-field triplet 
sub-levels of p x = \,y y = § , p z = [35] . Within the min- 
imal model discussed above, this would lead to a high re- 
action yield anisotropy, two-thirds that of the maximum 
possible (see Appendix). 

The use of spin-polarized triplets should open new 
channels for the design of bio-inspired molecular devices 
for sensing the direction of weak magnetic fields. 

We thank DARPA (QuBE: N66001-10-1-4061) and the 
EPSRC for financial support. 

Appendix 

Basis states. The spin dynamics of radical pairs may 
usefully be described in terms of two distinct sets of basis 
states. In both, singlet and triplet states are eigenstates 
of the total electron spin operator, S: 

(S|5|S)=0 (9) 
(Ti\S\Ti)= 72 (i = 0, ±1 or x,y,z) 

The triplet basis states are either the eigenstates of S z , 
the component of S along the z-axis: 

(T m \S z \T m )=m (m = 0,±l) (10) 

or are defined in terms of the three cartesian components 
of S: 

(T q \S q \T q )=Q (q = x,y,z) (11) 
(T x \S y \T z ) — i (and cyclic permutations of x, y, z) 

The relations between the two are: 
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|S) and |T m ) can also be written: 

IS) = i [K/? 2 > - \p ia2 }} (13) 
|T+i) = \aia 2 ) 

|To) = ^ + |A«2>] 

|T_0 = |/3i/? 2 > 

where \atj) and |/3j) are defined by: 

(a^Jaj) = +1 (14) 
<J3 s \8j,,\l3 s ) = -\ (j = l,2) 
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The axis system in which S x 



S y and S z 



and Si z are 



defined may be chosen to be the 'laboratory frame', in 
which the z-axis is commonly the direction of the applied 
magnetic field, or a 'molecular frame', which could, for 
example, be the principal axis system of one of the hy- 
pcrfine (HF) tensors. The triplet state \T q ) (q = x,y,z), 
Eq. (Ill, is spin-polarized in the q = principal plane 
within the molecule [19]. 

Initial state. The initial state of the radical pair spin 
system, p(0), is written as the direct product of the initial 
density operator for the two electron spins, po, and iden- 
tity operators for each of the nuclear spins (i = 1, 2, • • • ) 
to which the electrons are coupled: 



P(0) 




(15) 



(M is the total dimension of the nuclear spin-space). It 
is assumed that the formation of the radical pair is not 
nuclear spin-dependent. In the absence of chemical re- 
activity and spin-decoherence, the probability that the 
radical pair is in a singlet state at time t is given by the 
expectation value of the singlet projection operator, P s : 



P s )(t) = Tr p(t)P 



(16) 



Tr e 



-iHt 



p(0)e 



-imps 



where H is the time-independent spin Hamiltonian and 



P s = {|S)(S|}0 




(17) 



Minimal model. The principles of radical pair magne- 
toreception can be discussed using a simple model com- 
prising two electron spins and one spin-V2 nucleus (e.g. 
1 H) with an axially anisotropic HF coupling to one of the 
electron spins: 

HhR — ^ A qq Sl tq i q (18) 

x,y,z 

Si^ q I q + aa (§1 



q=x,y,z 

a 



a is the isotropic part of the HF interaction (expressed as 
an angular frequency) and a is a dimensionless axiality 
parameter [11]. Defined in this way, the HF interaction 
has cylindrical symmetry around the molecular z-axis. 

Two cases are considered specifically: a = (isotropic 
HF interaction) and a = — 1 (the HF interaction that 
results in the largest anisotropy in the reaction yield of 
this 3-spin system [T5]). 

The electron Zeeman interaction is included by means 
of the spin Hamiltonian: 



,7=1,2 



Sj, x sin 9 



(19) 
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FIG. 3. Reaction yield anisotropy, A$s, of the minimal 
radical pair model with a non-coherent initial state, po = ji. 
A<J>s is shown as a function of the rate constants fcs and k^. 
a = 1.0 mT, a = -1, w = 50 pT. 



in which u) is the strength of the applied magnetic field 
(expressed as an angular frequency) and specifies its 
direction with respect to the symmetry axis (z) of the 
HF tensor. It is assumed that the g-tensors of the two 
radicals are identical and isotropic, and that the nuclear 
Zeeman interactions are negligible. Both are excellent 
approximations for organic radicals subject to the weak 
magnetic fields of interest here. 

To account for the chemical reactivity of the radical 
pair within the minimal model, we use the 'exponential 
model' |23) in which singlet and triplet states react spin- 
selectively with the same first-order rate constant, fcg = 
fcx = k, to form distinct products. Although unlikely 
to be strictly valid for any real magnetoreceptor, this 
approximation simplifies the algebra without distorting 
the underlying physics. The yield of the chemical product 
formed via the singlet pathway, and its anisotropy, are 
calculated as l23l: 



P b )(t)e- kt dt 



A$s = max ($g) - mm ($s) 



(20) 
(21) 



so that < 4>g < 1. The corresponding yield for the 
triplet reaction channel is simply 1 — $g. $s is referred 
to as the reaction yield and A$g as the reaction yield 
anisotropy. The variation of $s with the orientation of 
the radical pair with respect to an external magnetic field 
forms the basis of the compass mechanism [2]. 

When fcg ^ fcx, the calculation of $g is performed in 
Liouville space: 



$ s = k s (P s \L- 1 \p(0) 



L = i [H (gills - Is® H 

+^k s (p s ® i 8 + i, 

+ \k T (P T ® i 8 + i 8 <8> P T ) 



(22) 
(23) 



P l 



G 



0.3 - 
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FIG. 4. Magnetic field effect on the reaction yield of the 
minimal radical pair model, k — 10~ 3 a. (a) Anisotropic HF 
interaction (a — —1). S (solid lines) and T (dashed lines) 
denote the initial radical pair states po = po(S) (p = 1) and 
po = Po(T) (p = 0), respectively, in Eq. J3|. The pertur- 
bation theory result in Eq. Q is valid in the shaded region 
where \a\ S> u 3> k. The dependence of $s on the strength of 
the applied magnetic field {oj/a) is shown for various angles 
between the symmetry axis of the HF tensor and the mag- 
netic field vector. The sharp features near log 10 (w/a) = 0.3 
arise from level anti-crossings (b) Isotropic HF inter- 

action (a — 0). S (solid line) and T (dashed lines) denote 
po = |S)(S| (v = 1) and po = |T z )(T z j (rj = 0), respectively, 
in Eq. |7|). The perturbation theory result in Eq. Q is valid 
in the shaded region where \a\ cj ^> k. The dependence 
of $s on the strength of the applied magnetic field (oj/a) is 
shown for various angles between the triplet alignment axis 
and the magnetic field vector. 



where is is the identity operator in the 8-dimensional 
spin-space and P T = Is — P S - 

Fig. [3] shows simulations for the minimal radical pair 
with an anisotropic HF coupling and an initial state: 



po = IMS) + fpo(T) = \i 



(24) 



A$g is non-zero except when kg = k^. The radical pair 
can exhibit magnetic compass properties even when its 



initial electron spin state is neither entangled nor coher- 
ent. The coherence arises during the spin evolution as 
a result of the differential reactivity of the singlet and 
triplet states. 

Perturbation theory. To obtain estimates of the max- 
imum possible magnetic responses within the minimal 
model, we use a perturbative approach [141 24J, appropri- 
ate for weak applied magnetic fields and long-lived radical 
pairs. This approximation is valid when \a\ uj k. 
These conditions are not unrealistic: HF interactions are 
of the order of 10 8 rad s _1 (sa 500 pT), the geomagnetic 
field is roughly 10 7 rad s^ 1 (« 50 f/T), and plausible 
values of the rate constant k are 10 5 — 10 6 s _1 [5]. 

Eqs ^ and Q were verified by the exact numerical 
simulations, the results of which are shown in Fig. [4j 

Multinuclear radical pairs. In the general case, the 
HF component of the spin Hamiltonian has the form: 



E E [ g j • A 3k ■ h 



(25) 



3=1.2 k 

= E E \_ a i^3 ' Ik + • T jk ■ ife 

3=1.2 k 

where ajk, Tjk and Aj^ are, respectively, the isotropic 
HF coupling constant, the anisotropic HF tensor and the 
total HF tensor for nucleus k coupled to the electron in 
radical j. 

The Zeeman term is: 



Zccman 



CO 



\Sj jX sin 9 cos ( 



(26) 



3=1,2 



5, 



where 9 and cj) specify the direction of the field in the 
molecular axis system. 

Fig. [5] shows, for completeness, versions of Fig. [T] in 
which the three initial triplet states are compared with 
|S)(S|. 

Concurrence. To quantify the entanglement of the 
various initial electron spin states po, we use the 'con- 
currence' proposed by Wootters for a two-qubit density 
operator [25] : 



C(p ) = max {0, A1 - A 2 - A 3 - A 4 } 



(27) 



where the Ai are the non-negative real square roots of the 
eigenvalues, in decreasing order, of: 



PO (°y (g) CFy) /So (fry <g> fry) 



(28) 



in which fry is twice the S y operator for a single electron 
spin, and /5g is the complex conjugate of po- 

General initial conditions. Within the minimal 
model, the most general initial state consistent with both 
Eqs Q and ^ is: 



p = £|S)(S| + (l-£ 



E 

q=x,y,z 



Pq\ T q)( T q 



(29) 
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(a) 0.15 
0.10 
0.05 

0.00 

(b) 0.15 

0.10 
0.05 

0.00 

(c) 0.15 

0.10 
0.05 
0.00 



TABLE I. HF data for N5 and H5 in FADH*. 






" 2 \og w 5 - 1 



FIG. 5. Reaction yield anisotropy, A$s, calculated for a 
radical pair in which one radical contains a 1 H nucleus (spin- 
y a ) and a 14 N nucleus (spin-1). k = 10 6 s" 1 and uj = 50 /xT. 
The HF coupling parameters (in mT) are as in the caption 
for Fig. [l] po = |S)(S| (black). p = \T*)(T X \ (a, green), 
= |T y )TT y | 



Po 



(b, green), p = |T Z )(T Z | (c, green) 



with < £ < 1, p x ,PyiPz > and p x +p y +p z = 1. Some 
such states are entangled and some are not. Almost all 
are anisotropic and may lead to compass behaviour even 
when the spin-Hamiltonian is isotropic. For the minimal 
model, with an isotropic HF coupling and \a\ 3> uj 3> k, 
the initial state in Eq. (29) leads to: 



$ s = I(2 £ +l) + i(l- £ ) 

2 a „„„2 



(30) 



cos^ (f> + py sin 2 9 sin 2 < 



(Px sin 

A$ S = z(l - e)(Pmax - Pmin) 



" Pz cos 2 6) 



Nucleus 


a 1 mT 


Tjj 1 mT 


Principal axes 




N5 


0.393 


-0.498 


0.4380 


0.8655 - 


-0.2432 






-0.492 


0.8981 


-0.4097 


0.1595 






0.990 


-0.0384 


0.2883 


0.9568 


H5 


-0.769 


-0.616 


0.9819 


0.1883 - 


-0.0203 






-0.168 


-0.0348 


0.2850 


0.9579 






0.784 


-0.1861 


0.9398 - 


-0.2864 




where p max = max{p x , p v , p z } and p min = 
miii{p x , p y ,p z }. The maximum anisotropy is ob- 
tained when e = 0, p max = 1 and p mul = (giving 
A$ s = i). 



The concurrence C(po) of the state in Eq. (29) is 

2e — 1 when e > | 
2p max (l - e) - 1 when e < | 

and p max > 
otherwise 



2(l-e) 



(32) 
(33) 

(34) 



Laboratory-frame polarization Although not relevant 
for a geomagnetic compass sensor, radical pairs can be 
created with large laboratory-frame polarizations, i.e. 
with unequal populations of the triplet eigenstates in a 
strong magnetic field (|T m ) (m = 0, ±1), as defined by 



(31) 



Eq. (10), with the z-axis being the direction of a strong 
applied magnetic field). For example, the Triplet Mecha- 
nism of Chemically Induced Dynamic Electron Polariza- 
tion can result in large polarizations for radical pairs pro- 
duced by triplet states formed by anisotropic intersystem 
crossing E3] ■ Another example, which also requires 
the electron spins to be quantized by strong electron Zee- 
man interactions, is seen in the EPR spectra of spin- 
correlated radical pairs with non-zero electron-electron 
exchange and/or dipolar interactions [5HI55] , 

Intersystem crossing. The implications of anisotropic 
intersystem crossing for the behaviour and properties of 
radical pairs is discussed in detail by Steiner and Ulrich 
[T7] (pp. 109-112). In most cases, as here, intersystem 
crossing is assumed to be independent of HF-couplcd 
nuclear spins. However, Kothe et al. [22 ], in an ele- 
gant study of quantum oscillations in an organic triplet 
state, have shown that the nuclei are in fact involved and 
that the appropriate molecular-frame triplet basis states 
are eigenstates of the combined zero-field and hyperfine 
Hamiltonians. We do not consider this possibility here. 
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FADH* radical. The HF interactions for FADH* used 
to calculate the reaction yield anisotropics shown in 
Fig. [2] are based on the data in Table [I] [28] . 
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